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\Q ' We prove theorems characterizing the minimizers in a model for condensation based on the 

Cahn-Hilliard free energy functional. In particular we exactly determine the critical density for 
\ droplet formation. 
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Let f2 be the d-dimensional square torus with volume L d . Consider the free energy functional J-(m) 
defined by 

F(rn) = — I |Vm| 2 dx + / F(m)dx , (1.1) 
2 Jn Jn 

where 9 is a parameter with the units of distance, and 

F(m) = ^(m 2 - l) 2 . 

This is a "double well" potential with minima at m = ±1. The two minima are the two "phases" 
of the system. The function mix) is an "order parameter field", representing a summary of the 
microscopic information about the underlying system locally at the point x that is necessary to 
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1 Introduction 

1.1 The variational problem 
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compute a density for the Helmholz free energy. The particular free energy density considered here 
is phenomenological; it does not arise from any particular microscopic model, but it is a simple car- 
icature of the Helmholz free energy densities that do arise from scaling limits of actual microscopic 
systems with phase transitions. See 0] for discussion on the relation with the microscopic models. 
The free energy functional is frequently called the Cahn-Hilliard free energy functional, but 
it has a fairly ancient history. It was already discussed by van der Walls JS] in the nineteenth 
century. 

For any number n with — 1 < n < 1, define the minimal free energy function /^(re) by 



/i(n) = inf | ^{m) : j-^ J rn(x)dx = n \ . (1-2) 



In what follows we shall work in units with 



This makes L and x dimensionless. That is, we are implicitly introducing dimensionless coordinates 
x, = x/0 and a dimensionless scale parameter L = L/9. However, setting = 1, we can drop the 
tildes, and work directly in dimensionless coordinates. 

If we fix a value of n with — 1 < n < 1, the minimization problem (|1.2|) simplifies as L becomes 
very large: In the bulk of £1 it must be the case that m(x) ~ ±1 to high accuracy. Consider any 
order parameter field m(x) that takes on only the two pure phase values ±1. Let V+ be the volume 
of the region in which m(x) = 1. While such an order parameter field would be discontinuous, and 
would therefore yield an infinite free energy, The quantity V + is quite relevant to (|1.2|) for large L. 
From the constraint m(x)dx = nL d , V+ — (L d — V+) = nL d so that V+ and n are related by 

v+ = n+± L d and n = _i + 2 Y+, (i.3) 
2 L 

Any order parameter field m{x) that is a minimizer for (|1.2j) will be continuous, and in fact, 
C°° as a consequence of the Euler-Lagrange equation that it must satisfy. Define the interface V 
between the two phases by 

r = { x : m(x) = } . 

Let |r| denote / H d ~ 1 (T), the d — 1 dimensional Hausdorff measure of T. Because of the gradient 
terms in the free energy functional, a minimizer m cannot make a sudden transition for +1 to —1 in 
crossing T, and so one would expect to pay a price for making such a transition that is proportional 
to |r|. That is, we would expect that 

f L (n) « 5|r| (1.4) 

for some proportionality constant S, which is called the surface tension. For the physics behind 
this terminology and this approximation, see |15j . 

Accepting this for the moment, if the volume of the region in which m(x) > is V, we would 
expect r to be a surface of minimal area bounding a region of volume V. The solution of the 
isoperimetric problem on the torus Q, depends on V/L d . Let denote the surface area of the unit 
sphere in W 1 . Since the volume of the corresponding unit ball in R rf is o~d/d 7 a volume of size V 
can be enclosed either by a sphere of radius (V '/ (o~d/ 'd)) l l d ', or in circular cylinder of radius s and 
length L. In d = 2, this is a "strip" between two parallel lines with length L each. Because we 
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are in the torus, the caps of the cylinder are not part of the boundary. Thus the surface of such a 
cylinder is <Jd-is d ~ 2 L, while the volume is {u d —i /d)S d ~ l L. Equating these to the surface area and 
volume of a sphere of radius r, a d r d ~ l and (ad/d)r d respectively, we see that both are equal, when 
r = r c where 

r c =((d-l)/d) d - 2 (a d ^/a d )L. (1.5) 

It pays to take advantage of the periodic boundary condition by reaching out to the boundary only 
for very large volume: V must exceed (a d /d)r d , which is proportional to L d , in order for this to be 
advantageous. For smaller volumes, the solution of the isoperimetric inequality is still given by a 
sphere. 

In what follows, we shall be concerned with values of n sufficiently close to —1, so that, with 
V = V+, as given by ()1.3j) . the minimal bounding surface is a sphere. 

The equimolar radius, tq, is defined to be the radius of the sphere whose enclosed volume is V+, 
where V+ is given in terms of n by (|1.3|) . Since (o~ d /d)r d = V+, we have from Q1.3|) that 

(±. ,.\"'. /v+^" d 



ro ^_(„ +1)j L . (1.6) 

Here we shall only consider values of n for which a d r^ 1 < 2L d ~ 1 , or, in other words, 

ro < r c (1.7) 

with r c given by fi.5|) . 

1.2 Two simple trial functions 

Under the condition (|1.7j) on n and L, there are two natural trial functions to consider for the 
variational problem (|1.2|) . 

The first of these is the equimolar droplet trial function 

m emd {x) = m (\x\ - r (n)) , (1.8) 

where we are representing as the centered cube in M. d with side length L and periodic boundary 
conditions, and mo{z) is some function such that linx^-i-oo 7710(2) = =Fl, with the transition from 
+ 1 to —1 being made in such a way as to minimize the cost in free energy. In fact, we require that 
the limits lim z _ i ± OQ mo(z) = =pl are achieved at finite values of z such that (|1.8|) does indeed define 
a smooth function on f2. 

The second of these is the uniform trial function 

m uni (x) = n , (1.9) 

corresponding to a "supersaturated" state with the order parameter strictly between the minimizing 
values. 

Accepting the validity of the approximation l|1.4jl . we have J'(m em( i) « Sa^r^ 1 , and therefore, 
from JHJ) 

^(m cmd ) « Sa d i^—(n + l)j L d ~ l = Sa d i^—^j , (1.10) 
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one easily computes that 

■n-uni) = \(n* - l?L d = ^ (l - • (1.11) 

Which of these trial functions provides a better description of the minimizers in ()1.2|) ? That 
depends on n, or what is the same, on the ratio Vj r /L d . In fact, there are now two obvious scaling 
regimes to consider: We can take L to infinity while keeping either n or V+ constant. For this 
reason, we have expressed .F(m un i) and J r (m em d) in terms of both n and V+. 

If one holds n constant, and takes L — > oo, then m emc j does much better than m un i. On the 
other hand, if one holds V+ constant as L tends to infinity, we see from and 1)1.10(1 that m un i 

does much better than m em( j for large L: This suggests that a droplet of the +1 phase will always 
"evaporate" into the surrounding —1 phase if the ambient volume |f2| is sufficiently large compared 
to V+. 



1.3 The critical scaling regime 

The situtation is much more interesting if one considers /z,(ra) with n tending towards —1 at the 
same time that L tends to infinity: we seek the smallest value of n(L) for which droplets are stable 
in a box of volume L d , and seek also to determine the structure of such critical minimizing droplets. 

In this sort of scaling regime, V+/L d will be very small, and we can express J-(m un i) in more 
physically meaningful terms as follows: Define the compressibility \ by 

X = ■ (1-12) 



Then since F(-l) = F'(-l) = 0, 

F(n) = F(-l + y + /L d )«-(^ 



1 fV^ 2 



2x \L d 

this gives us the approximation 

1 vl 

Hmum) ~ • (1-13) 

Of course, in our problem, x = 1/2- But introducing the compressibility highlights a competition 
between surface and bulk terms in minimizing the free energy. 
When Vl +l,d x L d , 

( v+ y- i/d i vi 

F(m em d) ~ Sa d [ — — and F(m UJ> i) « — — r 

\a d /dj 2xL d 

are comparable. For this reason, we refer to vl^ 1 x L d as the critical scaling regime. In terms 
of n and the equimolar radius ro, the critical scaling regime is characterized by 

n + 1 x L- d,{d+1) or, equivalents r X L d/{d+1) . (1.14) 



What should one expect for the minimizing free energy in the critical scaling regime, and will the 
minimizers be given by some sort of droplet, or not? 
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In a recent and incisive investigation of droplet formation in 2 dimensional Ising model 0], 
Biskup, Chayes and Kotecky proposed that to answer this question, one should introduce a volume 
fraction 77, and put 77V+ into the drop, and (1 — 77) V+ into the uniform background. They then 
constructed a phenomenological thermodynamic free energy function <&(r]) which is the sum of the 
surface tension term and the uniform background term: 

Here, < 77 < 1, and the suggestion in [I] is that in great generality, one can resolve a competition 
between surface and bulk energy effects by choosing 77 to minimize <3?. Defining C(n) by 



2 X S \dj I L d J d X S \d 
and I To I by |To| = <7d r o _1 > the quantity in (|1.15|) can be written as 

d>(r ? ) = 5|r |(r / 1 - 1 ^ + C (n)(l-r / ) 2 ) . 

Notice that 

*'")-*"" = ^ + C(nh - 2C ( B )), 

b\L 0| 

By the arithmetic-geometric mean, 

d/(d+l) 

> 



(1 1 -1 \ a/(a+ij 



Ql/(d+l)_ d + 1 



dd/(d+l) 



(1.17) 



Therefore, a minimum occurs at 77 > if and only if 

Let C* be the value of C that gives equality in this last inequality. One finds, as in @], 

„ 1 fd+i\< d+1 >' 2 , , 

ft = 3 (— ) • (118) 

Moreover, with C = C*, there is equality in the application made above of the arithmetic geometric 
mean inequality if and only if r]~ l l d /d = C*. Therefore, define 77* by 77* = (dC+)~ d . One finds 

— J ■ 

The heuristic argument of [I] suggests that when <& is minimized at 77 = 0, one puts all of the 
mass into the uniform supersaturated state, and there is no droplet. This is the case if C(n) < C+. 
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On the other hand, if is minimized at a strictly positive value of 77, then a strictly positive fraction 
of the mass should go into a droplet. This is the case if C(n) > C*. Moreover, it is easy to see that 
for all C(n) > C*, the minimizing value rj c of rj satisfies r\ c > 77*. As emphasized in 4 , this suggests 
that there are never drops containing a volume fraction less than 77*. That is, at least according to 
this heuristic analysis, there are never droplets whose volume is less than 

^{a d /d)r d . 

The validity of this was rigorously established for the 2 dimensional Ising model in [5]. We show 
here that the same heuristic analysis is correct for the minimization problem (|1.2j) concerning the 



f L 

Cahn-Hilliard free energy function T '. The first result concerns the value of the ratio — — -(n) for 

l r o| 

n = — 1 + K L~ d /( d+1 ) , for any K > 0, as L tends to infinity. 
1.1 THEOREM. For all K > 0, 

lim (-1 + KL- d ^ d+1 ^ = q mf i S (V" 1 ^ + D(K)(1 - r?) 2 ) (1.20) 



where 

D{K) = C (-1 + KL- d ^ d+1 ^ 



a d ^-i/d/K^^ d 



d x S \dJ V 2 



and S = 2 3 / 2 /3. 

Furthermore, let K± be defined by 



d J \ 2 J 

Then for all K < K*, and all L sufficiently large, the infimum in hl.2(J\) is a minimum attained 
uniquely at r] = 0, while for all K > K+, and all L sufficiently large, the infimum in hl.2(J\) is a 
minimum attained uniquely at 7] = r] c where r] c > 77*. 

To prove Theorem 11.11 we prove precise upper and lower bounds on /l(^) for values of n in 
the critical scaling regime, and from these bounds deduce (|1.2U|) . The remaining statements in 
the theorem then follow from the discussion just above concerning the minimization of ^(f?). For 
example, note that is obtained by solving D{K) = C* for K. The upper and lower bounds on 
fh{n) will be presented and proved in sections 2 and 3 respectively. We conclude section 3 with 
the proof of Theorem 11.11 

The theorem suggests that the curve n{L) = — 1 + K^L^ d ^ d+l ^ is critical for droplet formation, 
so that for large L and densities n significantly below this level, the minimizers will be uniform, 
while for large L and densities n significantly above this level, the minimizers will correspond to 
droplets of a reduced radius 7]l^ d rQ. The following theorems bear this out. 

1.2 THEOREM. For all K < K+ and L sufficiently large, when 

-1 < n < -l + KL- d ^ d+ V , 



the unique minimizer for hl.lfy is the uniform order parameter field m(x) = n. 
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Before stating the result concerning droplet minimizers, we must make this notion precise. To 
facilitate this, regard Q as the centered cube in M. d with side length L and periodic boundary 
conditions. 

For given r\ and n, and hence for given rj and tq, define a sharp interface reduced radius droplet 
order parameter field m^ n (x) by 



1 if \x\ < Tj^'^TQ 

-1 if \x\ > r/^ro 



1.3 THEOREM. For all K > K±, e > 0, and L sufficiently large, when 

-l + KL-V^ <n< -l + L- 1 ' 2 , 
any minimizer m for il.HAl is such that, after a possible translation on the torus Q, 



1 



m(x) - mj c n (x) 



4 

dx < e 



where r\ c is the minimizing value of rj in \1.2u\l . 

This theorem says that for large L, the set on which m and m|L n differ by an appreciable 
amount is small compared to (o^/d)rg. In particular, on the ball where m^ Cin = 1, m must be very 
close to 1 on all but a negligibly small percentage of the volume of that ball. Likewise, on the set 
where m\ Ctn = — 1, m must be very close to —1 on all but a set whose measure is a negligibly small 
percentage of the volume of the ball. The role of the fourth power is to make this small difference 
even smaller so that it is not overwhelmed by the large volume of the region external to the ball. 
Any power larger than 3 would work just as well in our argument. 

In this sense, m "looks like" m' Cjn for large L, and thus describes a droplet of the radius 
predicted by the heuristic argument of [1]. 

Theorems II .21 and 11.31 are proved in Section 4. Finally, in Section 5, we give an explanation for 
the remarkable efficacy of the simple trial function used in Section 2. In particular, we see which 
features of the free energy functional (|1.1|) are responsible for this. The point is that for the free 
energy functionals coming from models with a non local interaction, such as the ones considered 
in [S], as well as in |5] and these features are not present. However, the analysis in Section 5 
leads to a method for constructing trial functions of high accuracy that does apply to such cases, 
as well as to (jl.lj) . 
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2 The upper bound 

2.1 The interpolating family of trial functions 

For < rj < 1, let r„ = r/ 1//rf ro be the radius of a ball whose volume is rj times the volume of a ball 
with the equimolar radius, r$. The arguments of Biskup, Chayes and Kotecky suggest that one 
should use as a trial function a function of the form 



m n dx{x) = m (\x\ - r v ) + a(rj) , (2.1) 

where mo is a transition profile that very nearly minimizes the cost in free energy of making 
the transition from m = +1 to m = —1, and a(r]) is a constant determined by the constraint 
Ja m 'ndr(x)dx = n\fl\. As in Section 1, we are taking Q to be the centered cube in M. d with side 
length L. 

As r\ varies in the interval < r\ < 1, the family of "fractional droplet" trial functions defined 
in (|2.1|) interpolates between m un i, for r\ = and m em( j, for rj = 1. Of course, it remains to choose 
m . 



2.2 Planar surface tension and the choice of m 

The natural choice for mo is given by considering the problem of minimizing the cost per unit area 
in free energy of an infinite planar interface between the +1 and i phase. Denote this quantity by 
5; it will turn out to be the same constant S that appears in (|1.4|) . That is, 

5' = inf{ / ( -\m'(z)\ 2 + F(m(z)) ) dz : lim m(z) = ^1 
I Jm V 2 / z-»±oc 

Let m denote minimizer for this variational problem with m(0) = 0. The Euler-Lagrange 
equation satisfied by m is fh"(z) = F'(fh(z)). Multiplying both sides by fh'(x), and integrating 
from — oo to z, we obtain 

(fh'(z)) 2 = 2F(m(z)) , (2.2) 

since lim^-oo m'(z) = lim^-oo F(m(z)) = 0. 

One now easily deduces that fh(z) = — tanh(z/\/2), from which one could compute 5*. However, 
there is anther route that is more informative and useful in what follows: From (|2,2[) . we see that 



S = 2 F{fh{z))dz , 

J — oo 

and furthermore, 

w ,u, r° F(m(x)) , , 

F[m(x))dx = I — - — m (x)dx 

> J-oo m'[x) 



— oo 



F(m{x)) , 
m (x)dx 



'■man 



(2.3) 
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Thus, 



S = J ^2F(h)dh , 



(2.4) 



and hence S = 2 3//2 /3, which provides the numerical value for S that is quoted in Theorem 11.11 
However, in what follows, it is the integral formula, and not so much the numerical value, that 
turns out to matter. 

We now choose mo- We cannot simply choose mo = m since then mo(|x| — r v ) would not define 
a smooth, or even continuous, function on £1. However, only mild modification are required. Since 
we are interested in values of r with r = O (^L d ^ d+1 ^ , define 



m (z) 



fh{z) if \z\ < L( d - 1 )/( d+1 ) 
-sgn(z) if \z\ > 2L( d - 1 )/( d + 1 ) ' 



and smoothly interpolate in such a way that mo, like m, is odd. With any such interpolation, 
m o(M — r r)) defined a smooth function on and the difference between mo and m goes to zero 
exponentially fast as L tends to infinity. 

2.3 The determination of a(rj) 
The constraint equation is 



m ric j Lr (x)dx = nL d = 2(a c [/d)r d — L 



and hence 



a( V )L d = 2(a d /d)r d -L c 



mo(M — r^dx . 



(2.5) 



We require sharp estimates on the integral on the right. 



2.1 LEMMA. Define the constant M by 

M = I (sgn(z) - tanh(z/\/2))zdz . 
Jr 

For all L^- 1 )/^ 1 ) < r v < r , 

m (\x\ - r v )dx = [2(a d /d)r d V - L d ] - (d - l)M^- 2 ^ 2)/d + O (iT^) 



when d = 2 or d = 3. For higher dimension, the only difference is that the error term is 0(r d 4 ) 
Proof: Note that 



m (\x\ - r v )dx = 2{a d /d)r d - L d - I (sgn(|x| - r n ) + m (\x\ - r v ))dx 



Define 



I\ = I (sgn(|x| — r v )+mo(\x\ — r v ))dx and I2 = / (sgn 

J\x\<2r v ' J\x\>2r v 



(\x\—r v )+mo(\x\ -r v ))dx 
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We easily see that for all dimensions d, I2 = 0(e Ll/i ). Moreover, using polar coordinates, 



r 2r v 

h=°d I (sgn(s - r v ) + m (s - r v ))s d ~ 1 ds . 
Jo 



Introducing the new variable z = s — r„, we see that if we extend the integration in z over the 
whole real line, we only make an error of size 0(e~ Ll 4 ) at most, and so 



h = °dr d - X [ (sgn(z) + mo(z)) U - -) dz + 0(e" Ll/4 ) . 

Taking into account the fact that (sgn(z) + mo(z)) is odd and rapidly decaying, we see that for 
d = 2 or d = 3, 



/ (sgn(z) + mo(z)) (l - — J dz = - — - j (sgn(z) + m (z))zdz . 
Jr \ r vJ r v Jr 

In any dimension, this gives the leading order correction. This, together with the definition of 
mo(z) in terms of m(z) = — tanh(z/\/2), yields the result. □ 

Therefore, (|2.5jl together with Lemma 12 . 1 1 yield for d = 2 or d = 3 that 

a( V ) = 2(<7„/d)^(l -v) + (d- l)Ma d ^-r,^ d + 0(e~ Ll/A ) . (2.6) 

The only difference for higher dimensions d is that 0(e~ Ll/4 ) must be replaced by 0{r d ~^ / L d ). 
Notice that unless 77 = 1, the first explicit correction is already very small compared to the leading 
term; it is smaller by a factor of r^ 2 . In the critical regime, by (fmf) . 2 x L- 2d / d+1 \ Moreover, 
we see that in the critical scaling regime, except when 77 = 1, 



o 



(77) x . (2.7) 



2.4 Computation of ^(m^dr) 

With the trial function specified, we now compute ^{m^). 
2.2 LEMMA. In the critical scaling regime ro x L d ^ d+l \ 



„3d 

(d 2 -3d)/((i+l) 



^Kdr) < - 8(a d /d) d ^(l - r?) d + O \ LS a j , (2.8) 

w/iere £/ie /irsf term on £/ie n#fa is ^L( d2 - d )/( d+1 )^ and the second is O ^L( d2 ~ M )/( d+1 )^ . 

Notice that the leading term in the upper bound is exactly $(77), and that the next term is 
negative. 

Proof: To simplify the notation, we write too to denote mo(|a;| — r„) and a to denote a(rj) so that 
rn^dr = m o + a - Then 

F(m vdl ) = F(mo) + F'(TO )a + -F"(m Q )a 2 + V"(mo)a 3 + 7a 4 . 

2 4 
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We are required to produce a close upper bound on the integral of each of these terms over Q. 
We start with a f n F'(mo)dx. 

Note that F'(m) = m 3 — m, Since m^z) — mo(z) is an odd, rapidly decaying function of z, 
estimates just like the ones employed in the proof of Lemma (|2. 1|) show that 



d-l 



F'(m )dx = a d ri / (rffiz) - m(z)) 1 - - 

Jr V i 

Then, with the constant B defined by 

B = l (m 3 (z) — m(z))zdz , 
Jr 

we have for d = 2 or d = 3 that 



dz + 0(e 



f F'{m )&x = o d r*- 2 Bj 1 { - d - 2 V d + 0{e 
Jn 



(2.J 



and the same is true for d > 4 except that the error term must be replaced by 0(Vq 4 ). 
Next, F"(m ) = 3mg - 1 < 2 = 1/ X - Therefore 



Finally, F"'(m) = 6m, and so 



/J F » (mo)dx <_L L ,. 

/ fl-^ ?/ "( m o)dx = / modx 



(2.10) 



(2.11) 



and this integral has been computed in Lemma 12.11 
Now let / denote the integral 



I 



F'(m Q )a + lF'"(mo)a 3 + ]a 4 
o 4 



dx 



In the critical scaling regime, the dominant contribution to I comes from the F'" term, and is 
-L d a 3 >c L,( d2 - 2d )/( d + l ) . Each of the terms in the integrand contributes at the order ^(^-^^/(d+i) 
in the critical scaling regime, and we have 

I = -L d a 3 + [a^B^^a + 2(a d /d)r d V a 3 + a 4 L d /4\ + O (l(*-*»)/(*U)) . 
Using 1)2. 5 Jl . we can express this as 



+ 



-8(a d /df^ M (l-r ] f + 

q ^ L (d 2 -4d)/(d+l)^ 



Jld-2 Ad 

2d{o d /d) 2 B^- r r 1 2 - 2 l d + A(a d /d)^(l - V f(l + 3 V ) 



(2.12) 



where the first term on the right is proportional to , and the second is proportional 

tQ L (d 2 -3d)/(d+l)_ 
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Finally, we have to estimate L [|Vmo| 2 + F(mo)] dx. Once more, estimates just like the ones 
employed in the proof of Lemma 12.11 show that 

f [|Vm | 2 + F(m )] dx « a^ 1 [ [\m'(z)\ 2 + F(m(z))} (1 + z/r^dz 

where the errors are exponentially small in L 1 / 4 . But because mo is so close to to, this only 

differs from S\To\r] 1 ~ 1 / d by errors that are (D(r d ~~ 3 ). In the asymptotic scaling regime, r d ~ 3 x 
L {d?-M)/(d+l)_ 

Combining estimates, we have 



SlToK^ + C(n)(l - r/) 2 ) + O ^^)/(d+i)^ 



(2.13) 
□ 



3 The lower bound 

3.1 An A priori pointwise upper bound 

Standard compactness arguments show that the infimum in ()1.2|) is attained at a minimizer m{x) 
which satisfies the Euler-Lagrange equation 

—Am(x) + m 3 (x) — m(x) + [i = , (3-1) 

where \x is a Lagrange multiplier corresponding to the constraint in (J1.2j) . 

Our immediate goal is to prove an a priori pointwise upper bound on a minimizer m that is very 
close to 1 in the critical scaling regime. Such a bound can be obtained from the Euler-Lagrange 
equation and the maximum principle. 

Let 

^min and x max be such that for all x, 

m(x min ) < m(x) < m(x max ) . 

These exist since any solution of the Euler-Lagrange equation is continuous. 

We will now show that m(x mSuX ) cannot be too large. Define numbers A and v by 

1 + A = m(x max ) and - 1 + v = m(x min ) . (3.2) 

It will be convenient in the arguments leading to the proof to write n in the form 

n = -l + S. (3.3) 

Notice that in the critical scaling regime, 5 x L- d ^ d+1 \ Also, from JOJ) and (fOjl . 

S = 2^ = 2(a d /dA . (3.4) 
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3.1 LEMMA. For any solution of the Euler-Lagrange equation let A and v be given by 

\3. 6 J\) . Then v>\. Consequently, if m is any minimizer for (1), 

5 > A . (3.5) 

Proof: Evidently, Am(x max ) < 0, and so from (|3~T|> and (f3~2|) . (1 + A) 3 - (1 + A) + fi < 0, or 

2A + 3A 2 + A 3 < -n 
In the same way, from ()3.1() and (|3,2j) we have 

2f - 3z^ 2 + v 3 + n > -fi . 

Therefore, 

2u - 3u 2 + v 3 > 2A + 3A 2 + A 3 . 

It evidently follows that 

2u + v 3 > 2A + A 3 , 

and since /(x) = 2x + x 3 is monotone increasing, it follows that v > A. 

Next, since the average value of any function is no less than its minimum, it follows that 
n > — 1 + v, and by (4), this means 6 > v. Combining estimates, we have (|3.5|) . □ 



3.2 An A priori lower bound on m(x max ) 

We next show that any for non constant minimizer m, it cannot be that m(x max ) is much smaller 
than 1. For this purpose, define w by w(x) = m(x) — n. For m satisfying the constraint in (|1.2j) . 

w(x)dx = . (3-6) 

Clearly, J n |Vm| 2 d3; = J Q \Vw\ 2 dx, and 

/ \{m 2 - l) 2 dx = L 2 (n 2 - l) 2 + / - ( ^(3n 2 - l)w; 2 + nw 3 + ]w 4 ) dx , 

since the terms linear in w drop out due to (5b), and 

nw H — w = -w (w + 2n) — n w . 
4 4 

Hence, if we define the functional Q by 

Q(w) = - [ |Vtt;| 2 dx+ / G(u;)dx , 
2 Jn Jn 



where 



(3.7) 



G (w) = 1 ^-^w 2 + ^w 2 (w + 2n) 2 , (3.8) 



we have 

F(m) = F{n)+Q{w) (3.9) 
whenever m satisfies the constraint in (|1.2I) . 



CCELM May 9, 2005 



14 



3.2 LEMMA. Let m be any minimizer for \1. 6 J\) . and suppose that m is not constant. Then 

m(x max ) >l-5- 2^/5^1-5/2 . (3.10) 

Proof: Notice that 

n 2 — 1 1 . , 9 

— — + ~(w + 2n) 2 < 
2 4 

if and only if 

z_ < w < z + 

where 

z± = -2n ± y/2 - 2n 2 = 2 - 25 ± 2\/*\/ 1 ~ V 2 ■ 

Since m is not constant f n \Vw\ 2 dx > 0. Thus Q(w) > unless w(x max ) > 2 — 25 — — 5/2. 

If m is a minimizer, Q{w) > is impossible, on account of (|3.9jl . Since m(x max ) = n + w(x max ), we 
have the estimate. □ 



3.3 A partition of Vt 

We now partition 0, in to three pieces, one of which will contribute a surface tension term to the free 
energy, another of which will contribute a compressibility term, and another that will be negligible. 

Returning to the original dependent variable m, we fix a number k > to be determined below. 
However, to fix our ideas for the time being, suppose that k = 0(5 1 ^ 3 ). Define numbers h + and 
h- by 

h+ = 1 — k and h- = —1 + k . (3-H) 
Define the sets A, B and C by 

A = { x : h- < m(x) < h + } B = { x : m(x) < /i_ } , 

and 

C = { x : m(x) > /i+ } . 

If m is a non constant minimizer, and k = O^ 1 / 3 ), then for L large enough, C will be non empty 
by Lemma 3.2. Define a radius R by 



(o- d /d)R d = \C\, 



(3.12) 



where the right hand side denotes the measure of C. Evidently R is the radius of the ball with the 
same volume as C. 

We shall obtain a lower bound on fh{n) by separately estimating the integrals 



Ia 



-|Vm| 2 + F(m) 
2 v 1 



dx 



and 



Ib 



B 



^|Vm| 2 +F(m) 



dx . 



(3.13) 
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3.4 The surface tension contribution 

We now prove a lower bound on I a, which corresponds to the surface tension contribution to the free 
energy. The lower bound is obtained through use of the co-area formula pP, which expresses 
the volume element in $7 as 

dx = — — —-7da h dh , 

\\/m{x)\ 

where dah is the surface area along Fh, the level set {m(x) = h}, or, more properly put, the d — 1 
dimensional Huasdorf measure on this set. 

It is worth noting at this point that the rearrangement inequalities of the sort discussed in [H] 
apply in this case, and allow us to conclude that for any minimizer m, the level sets are symmetric 
monotone. If we translate so that the maximum of m is at 0, then this means that for any h, and 
any of the standard basis vectors ej, j = 1, . . . , d, the set of t for which m(tej) > h is a symmetric 
interval. In particular, is a rectifiable, simply connected curve. 

In applying the co-area formula, we shall gloss over certain standard technical issues. These 
are all explained, for example, in the discussion of the Faber-Krahn inequality in [7] or ^U], where 
the co-area formula is applied to another variational problem, namely the one for the fundamental 
eigenvalue for the Laplacian in a a domain in M. d . Those readers who are not familiar with the 
use of the co-area formula in proving inequalities such as the Faber-Krahn inequality may wish to 
consult the references cited above. For those who are, we proceed with the proof. 

3.3 LEMMA. Let m be any non constant minimizer for il.ty) . and suppose that k = OiS 1 / 3 ). 
Then for L large enough, Then 

I A > (o- d )R d - x (S - 2k) . (3.14) 

Proof: By the co-area formula, 

Ia = [ k+ [ (\\Vm{x)\ + t^%) da h dh . 
Jh- Jr h V 2 \Vm{x)\J 



By the arithmetic-geometric mean inequality. 

1 / F(h) 



and therefore, 



Ia> f + \T h W2F(h)dh 

Jh- 



where |r^| denote the one dimensional Haussdorf measure of T^. 

Note that encloses a region whose volume is at least \C\. By the isoperimetric inequality 
on the torus, the length of the boundary of such a region is at least a 1 J d {dV) l ~ 1 / d provided \C\ < 
(2 d /( d + 1 )x 2 )/(cfer^ _1 ), and is at least 2L d ~ x otherwise. ( If U is any domain in n with \U\ < 
(2 d /^ 1 )L 2 )/(do- d r 1 ), then the surface area of the boundary of U is no less than that of a ball of 
volume \U\, namely cr l J d (dV) l ~ l l d ). Otherwise, the lower bound on the perimeter is simply 2L d ~ 1 , 
and this is achieved by a "strip" of width \U\/L d ~ 1 in the torus.) Therefore, defining P by 

P = min{ a 1 J\d\C\) 1 - 1 ' d , 2L d ~ x } , 
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we have {T^l > P for all h- < h < h + . Hence 



I A >P y / 2F{h)dh 
Jh- 

By ((231), tnis y ields 

I a >p(s- J s/2F(h)dh - J y/2F(h)dhj 

Furthermore, 

^/2F{h) = (1 - h 2 )/V2 < 1 , 

and hence v / 2F(/i)d/i < k and ^2F(h)dh < k. This gives us Jyi > P(S — 2k). Now if 
|C| > L 2 /tt, this would imply i"^ > 2L(S" — k), which is much larger than J-(n). This is therefore 
impossible when m is a minimizer for ()1.2j) . and so with i? defined by ()3.12() . we have the bound 
<IT1 . 



3.5 The bulk contribution 

In this subsection, we prove a lower bound on the contribution to the free energy from B. For this 
purpose, we first require an upper bound on \A\ which shows that, for large L, \A\ is negligible 
compared to \C\. Ideally, one might hope that A is an annular region about C, and to obtain a 
"surface term" type bound for A, showing that is is bounded by a multiple of ?~q -1 x L d /( d + 1 ), 
However, as one can see from the proof of Lemma I3.3| even if C were spherical and A was an 
annulus about it, one would have to take the annulus to be fairly "thick" in order to capture most 
of 5 in the estimate (j3.14j) . Thus, the following simple estimate is rather sharp. 

3.4 LEMMA. Let m be any minimizer for \1. 6 J\) . Then 

jd r2 

\A\ < 2F(n)^ < 2- 2 L d . 

Proof: Since 

F(h+) = k 2 (1 - k/2) 2 = F(h_) , 
it is easy to see that uniformly on A, 

F(m(x)) > k 2 {1 - k/2) 2 (12) 

Therefore 

Ia > \A\k 2 {\- k/2) 2 . 
On the other hand, since m is a minimizer, 

I A < = F{n)L 2 . 



In the range of 5 being considered, (1 — k/2) 2 > 1/2. 



□ 
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Our next goal is a lower bound on Is- Notice that on (— oo, h~), F is strictly convex. In fact, 
F"(h) >3h1-l. Define the quantity X- by 

— = F"(h) = 3hl - I . 

X- 

Then, by Taylor's Theorem, and using the fact that F(— 1) = F'(—l) = 0, we have 

F(m(x)) > -—(m( x ) + l) 2 

everywhere on {m < h-}. 
Therefore, 

/ F(m(x))dx = \B\ ( - — - I F(m(x))dx ) 
Jb \\ b \Jb ) 

(m(x) + l) 2 dx 



> \B\ 

> \B\ 



1 ( 1 



2X- Vl^l 
1 ( 1 



'2x- \\B 



(m(x) + l)dx 



— — r ( J m(x)dx + \B\ 
-\ B \ \Jb 



b 



O 



We now need an upper bound and lower bounds on \B\ and J B m{x)dx. Note that 

m{x)dx = nL d — f m(x)dx — f m(x)dx . 
By Lemma 13. II and the definition of R, 

(1 - K)(a d /d)R d < I m{.r)d.r < (J + S)(rr,,/<I}R'> 

By Lemma 13, 41 



c 



Thus since k > 5, 



\A\ < h_\A\ < / m{x)dx < h + \A\ < \A\ . 
J A 

< \A\ + K{a d /d)R d 



< 2\A\ + K(a d /d)R d 



m(x)dx - (nL d - {<r d /d)R a 

Next, it is evident that \B\ = L d - (a d /d)R d - \ A\. Therefore 

(^j m{x)dx + \B\^j - (5L 2 - 2{a d /d)R d ) 

Hence, if we define e by 

e = 2(2\A\/L d + K(a d /d)(R/L) d )\S - 2{a d /d){R/ L) d \ , 

we have 

j m{x)dx + \B\^j > {5L d - 2{a d /d)R d f - eL 2d . 



(3.15) 
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By Lemma ITU \A\/L 2 = 0{5 2 /k 2 ), and for any mmmuzer, we must have R = 0(L d ^ d+1 \ since 
otherwise, if R were any larger, the contribution from the interface as estimated in Lemma fe. 31 would 
already exceed the free energy for the uniform trial function. Thus, (R/L) d = 0{L- d K d+l ">) = 0{$). 
Therefore, 

€ = O ( — g ~\~ k5 2 

With the choice k = <5 1//3 , this gives us 

e = 0{5 7/3 ) , (3.16) 

with the essential point being that this is negligible compared to 5 2 as L tends to infinity in the 
critical scaling regime. 

Finally, since \B\ < L d , this proves the following bound: 

3.5 LEMMA. Let m be any non constant minimizer for Then, with e given as above, 

L B >^-[(5-2(a d /d)(R/L) d ) 2 -e) . 



2X- 

It now follows from Lemmas 13.31 and 13.51 that for any non constant minimizer m, 

T{m) > Ia + Ib 

> a^- 1 (S - 2k) + ^-{8 - 2(a d /d)(R/L) d ) 2 - ^—e . 

2x- 2x- 



(3.17) 



It now remains to optimize this over R. 

Now introduce S- = (S — 2k) and r\ = R d frfi. Then we can rewrite this lower bound as 

T{m) > 5_|r | ( r } 1 - 1 ' d + _^-C(n)(l - rj) 2 ] - ^—e . (3.18) 



S-X- v " "J 2 X 

Proof of Theorem 11.11 As L — > oo in the critical scaling regime, 5_ — > S and X- ~* X- 
Moreover, from dS3SD , (j3~l) and the definition of e, L d e/\T \ -> as L -> oo. Thus, (I3~T%1) 
provides the lower bound needed to prove (|1.2(Jj) . The upper bound is provided by Lemma 12.21 
The remaining statements follow from the analysis of the minimization of the phenomenological 
free energy function (|1.15|) that was explained in the introduction. □ 



4 The structure of the minimizers 
4.1 The proof of Theorem 11.21 

Suppose that n = —1 + KL d l <yd+1 ^ where K < K*. We would like to conclude from (|3.18|) that any 
non constant trial function m has a higher free energy than the uniform trial function m(x) = n, 
at least for all sufficiently large L. 
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Recalling that S\To\C(n) = F(n), define r\ by 

V = sup{ V : 5_|r |(r ? 1 -^ + _^ C '(n)(l-r ? ) 2 )-^e<5|ro|C'(n)| 

As in the proof of Theorem ll.il for all L sufficiently large, S- is sufficiently close to S, and X- 
is sufficiently close to \ that 

for some C < C*. For C < C+, the unique minimizer of 

v ^l-i/d + C(l - r/) 2 

is 77 = 0. Therefore, since eL rf /|ro| — > as L —> oo, it follows that f\ — ► as L — » oo. 

Now, as in the previous section, for any non uniform minimizer m, there is a relation between 
r\ and the size of the level set \{m > 1 — k}\ given by w = (R/ro) 1 l d and \{m > 1 — k}\ = {o~d/d)R d . 
Here, as in the last section, k = with 5 given by (|.3..3|) . It follows from (|.S.18|) and the definition 
of f\ that for any non constant minimizer m, w < fj, and so \{m > 1 — k}\ is negligibly small 
compared with the volume of the equimolar ball; that is, (<Jd/d)rQ, when L is large. 

In other words, if n = — 1 + KL d <( d+1 ' where K < K+, and L is large, then any droplet in any 
minimizer must be extremely small. To prove Theorem 11.21 it therefore suffices to show that such 
extremely small drops are impossible in a minimizing order parameter field. We do this in the next 
lemma. 

4.1 LEMMA. For all K > 0, there is a constant Ck > depending only on K so that if 
n < — 1 + KL~ d /( d ~*~ 1 > and m is any non uniform minimizer for kl.Ofy . then 

\{m > 1-k}\ >C K r d . 

Moreover, Ck is uniformly strictly positive for all K in an interval around . 

Proof: We again work with the functional Q{w) which is defined in (|3,7j) and related to J-{m)—J-{n) 
by ([3.9)1 . Clearly, if F(m) < F(n), then the potential G(w), defined in (|3.8j) . must become negative. 
However, as seen in the proof of Lemma I3.2| G(w) < if and only if Z- < w < z + where 
z± = 2 - 25 ± 2VSs/l - 8/2. 

Moreover, G(w) > (n 2 — l)w; 2 /2 for all w. Since w = m — n and by Lemma 13.11 m < 1 + 5, 
while by definition n = —1 + 5, w < 2. Therefore, 

G(w(x)) > -45(1 - 5/2) 

for all x. 

Define the set C by C = {x : w(x) > z_ }, and define the number R by 

(a d /d)R d = \C\ . 

We have 



^\Vw\ 2 + G(w) 



dx> [ G(w)dx > -45(a d /d)R d . (4.1) 
Jc 
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Define the set A by A = {x : < w(x) < z_ } where the lower bound is arbitrary but 
convenient. The same argument used to prove Lemma 13.31 shows that 



^\Vw\ 2 + G(w) 



dx > SadR 1 - 1 , (4.2) 



i A L 

where S = /*" ^/2G{Kj&h. 

It now follows from ()3.9|) that J-(m) > J~{ri) unless 

A5(a d /d)R d > SadR 1 - 1 . 

However, J-(m) > J~{n) is impossible if m is a minimizer. Hence, R > S/(45). From (|3.3|) . 
5 = 2(a d /d)(r d+1 /L d ) — . Hence 

ft - Ld 

R -8(a d /d) r d Q + ir ° • 

The condition n < —1 + KL~ d ^ d+1 ^ yields the bound 

Thus, there is a constant Ck depending only on K so that R > CkTq. 

The final observation to make is w > z_ if and only if m > 1 — 5 = 2y/Sy/l — 5/2. Therefore, 
since k = (5 1 / 3 , C C C for all L large enough. □ 

Proof of Theorem 11.21 By Lemma 14.11 whenever m is a minimizer for (|1.2|) with n < — 1 + 
KL—d/id+i) an d j( < fag corresponding value of rj is bounded away from zero by a strictly 
positive quantity depending only on K. But by the remarks preceding Lemma 14.11 the rj value of 
any minimizer cannot exceed fj, which tends to zero as L increases. Hence, for L sufficiently large, 
there are no non constant minimizers with n < — 1 + KL~ d ^ d+1 ^ and K < K±. □ 

4.2 The proof of Theorem ITT31 

The essential tool here is quantitative version of the isoperimetric inequality. The classical version, 
valid in d = 2, is due to Bonnesen 0, who worked in the setting of convex geometry. For the 
form used here, see |T3]. To formulate the inequality, let U be a domain in M. 2 bounded by a 
simply connected rectifiable curve Suppose that p ut(^) is the imnmum of the radii of circles 
containing U, and Pi n (U) is the supremum of the radii of circles contained in U. p ut(U) is called 
the outradius of U, and Pm(U) is called the inradius of U. The inequality of Bonessen states that 

|r| 2 - Att\U\ > 7T( Pout (U) - p in (U)) 2 (4.3) 

where |T| denotes TL 1 ^), the 1 dimensional Haussdorf measure of T. That is, 

in ^L+je-mzem) . (4.4) 



\U\ ~ \ s/\U\ 
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We shall show how this inequality may be used to prove Theorem II .31 for d = 2, and then shall 
explain how a recent extension by Hall ^2] of Bonnesen's inequality to higher dimensions yields 
the result for d > 2. The essential points are clearest for d = 2, and so we begin with that case. 

We first have to justify the application of this inequality on the torus. The point is that by 
properties of the rearrangement employed in [SJ, we know the level sets of minimizers m must be 
connected. If they "wrap around" the torus, their perimeter has a length of at least 2L, which 
would give too large a surface contribution to the free energy 

To apply (|4.4|) . we return to the proof of Lemma 13.31 and take U = Uh = {m > h } for 
h- < h < h + . Then, using the notation of section 3, 

\U h+ \ = \C\ and \U h _\ = \C\ + \A\ . 

The essential point is that in the critical scaling regime, \A\ is negligibly small compared to \C\ for 
large L. This is the content of Lemma 13.41 Therefore, for all e > 0, if L is sufficiently large, 

\U h _\ < (l + ef\U h+ \ = (l + e ) 2 |C| . (4.5) 

Recall that in section 3 we have defined R by ~kR 2 = \C\. Now suppose that 

Pout(U h+ ) > (l + 2e)R . (4.6) 

Then since for all h < Uh + C Uh, 

Pout (£4) > (1 + 2e)i? for all h_<h<h + . (4.7) 
On the other hand, since for all h < h+, Uh C Uh_ , so that 

\U h \ < \U h _\ <7T(l + e) 2 R 2 , 

it follows that 

pin(U h ) < (1 + e)R for all h_<h<h + . (4.8) 
Combining (|4.4|) . (|4.7j) and (|4.8|) . and letting T^ denote the boundary of Uh, we have 

|r => v / vi : (2 + ^— J for all h_<h<h + . (4.9) 



Thus, under the hypothesis (|4.6|) 



|r A | > 2vr ( 1 + 2 + 2 J R for all h^<h<h + . (4.10) 

Going back to the proof of Lemma l3.3[ one see that effectively, the hypothesis ?? would introduce 
and extra factor of (1 + e/(2 + 2e)) into the surface tension S. This would increase the leading 
order contribution to the free energy over what what obtained in Section 2 with the trial function 
m^dr- Therefore, for all sufficiently large L, the hypothesis (??) is incompatible with m being a 
minimizer. 

The same conclusion can be obtained from Hall's theorem i n higher dimension, in essentially 
the same way. The version of Hall's theorem found in Theorem 4.1 of ^1] is particularly useful for 
this purpose. We summarize the discussion in a lemma, using the notation introduced above. 
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4.2 LEMMA. Let K > K± and n > -l + KL~ d ^ d+1 \ For any e > and all L sufficiently large, 
if m is any minimizer for \1. 6 J\) . then with B(xq,R) denoting the ball of radius R centered on xq, 
there is a point xo with 

\CAB(x ,R)\ 
\B(x ,R)\ ~ £ ' 

where A denotes the symmetric difference C U B(xq, R)\(C n B(xq, R)). 

Proof of Theorem 11.31 Let m be any minimizer. After translating, we may assume that 

\CAB(0,R)\ 



Also, since K > K+, we know that the value of R must be very close to nV d rQ. Therefore, for all 



\B(0,R)\ 

mow that the va 

L large enough, we have 



ICA^O^ro)! <2g 
15(0,^7-0)1 

Now, if \m{x) — mt, v (x)}\ > 2k, it must be that either x G CAB(Q,rjc d ro), or else x S A. 
Hence, the set 

{ x : \m{x) — > 2k } 

has a volume no greater than 

2e| J B(0,r )| + |,4| , 

and we recall that by Lemma 13.41 |j4| is negligibly small compared to r d for all sufficiently large L. 
Thus, we obtain 

|{ x : \m(x)-m{ v (x)}\>2K}\<3s\B(0,r )\ 

for all sufficiently large L. Of course we have the globally valid bound \m{x) — m n ^{x)}\ < 2. 
Therefore, 

/ \m(x) - ml (x)}\ 4 dx < 48e|B(0,r )| + 16K 4 L d . 
Jn 

Recall that k = 5 1 / 3 , and that in the critical scaling regime, 

6L d _ L d/(d+i) _ |s(0,r )| • 

Therefore, for any e > 



1 



\m(x) — mr n (x)}| 4 dx < e 



\B(0,r )\Jn f " ' ^ 
for all L sufficiently large. □ 



5 The construction of good trial functions 

5.1 A Chapman— Enskog— Hilbert expansion approach 

As we have seen, the simple trial function was sufficient to provide the upper bounds required 
here. In fact, it is quite likely that the upper bounds computed in Section 2 are accurate to at 



CCELM May 9, 2005 



23 



least the first two orders in powers of L. Our goal here is to present a systematic construction of 
high order trial functions. This will explain the remarkable efficacy of the simple prescription m v dr 
for , but will also suggest how one should construct trial functions of similar efficacy for other 
free energy functionals. To keep this section concise, we only treat the case d = 2. This is fully 
representative, except that the formulas are much simpler. 

For d = 2, in the critical scaling regime, ro x L 2 / 3 . We introduce the scaling parameter 
A = r-Q 1 . After the rescaling x' = Ax, so that U —* tl x = so that |S1 A | = L 2 A 2 = L 2 r^ 2 , the 
Euler-Lagrange equation (|3.1|) becomes 

-A 2 Am + F'(ra) + /i A = 0. (5.1) 

We will drop the prime on the new coordinate x for sake of simplicity. 

An approximate solution of order N to the Euler-Lagrange equation is a function mS N ^ such 
that 

-A 2 AmW + /( m W) + M W = 0(X N+1 ) (5.2) 
Since we are interested in solutions m to the Euler-Lagrange equation such that 

^ m(x)dx = n 



\n x \ 



we require that the constraint on the mass be approximately satisfied in the sense that the approx- 
imate solution m( N ) satisfies 

/ m (N \x)dx = n + 0(X N+1 ) . (5.3) 
I s * I Jn x 

Our aim here is to use an expansion method, based on the Chapman-Hilbert-Enskog expansion of 
kinetic theory to construct such approximate solution, and to use them as trial function for 1)1. 2j) . 
after adding a small constant so that the constraint is exactly satisfied. 

To do this, we first introduce local coordinates in a neighborhood of the curve r^, which 
will be determined in the course of the expansion. Let s denote an arc length parameter along 
pW. The starting point of the parameterization is immaterial. We denote by d(x, the signed 
distance of x from r^, with d(x,T^ N ^) > when x is in the interior of (i-e., the smaller of 
the two regions into which divides the tours). Define a "fast" variable z by z = d(x,T^ N ')/X. 
Then (s, z) give us a system of coordinates on a tubular neighborhood around 

To construct the approximate solutions, we make the following prescription, which has several 
parts. For any positive integer iV: 

(1) The interfacial curve will be a circle of radius r^ N ' to be determined at each order, 
essentially by the condition 1)5. 3 Jl . Note that, because of the rescaling the radius A N ' is actually 
measured in units ro and the in the original units the radius of the circle is r'^Vp 

(2) The chemical potential fj,( N ' has an expansion of the form 



(3) We then construct 



fiW = Ami + A 2 /i 2 + • • • + X N fi N . 



m (A0 _ m ( d(*,r {N \ , \ (5.1) 



\ N 

- + A " K + <Pn] 
/ n=l 
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where: (i) thq is the approximation to m introduced in Section 2, (ii) <pj will be a bounded 
continuous function that is nearly constant away from T^ N \ (Hi) hj is a function that has the form 



h; 



where on the right, hj(z) denotes a rapidly decaying function of the variable z (The notation is 
such that the symbol h n plays two roles, but this should cause no confusion.) 

By using the local coordinates (s, z) around the curve r^, we write the Laplacian as 

A 2 A f — ^— f -\ XKiN) 9 fl ^ 92 (55) 

and 

n=0 n=0 

Note that our simplifying assumption that r'^' is a circle amount to assume that K (7V) does not 
depend on s, thus dropping a few terms of higher order in A in (|5.5[) . Note also, for future reference 
that the area element in these coordinates is 

(1 + K ( - N hz)Xdzds (5.7) 

and recall that the surface tension, S, is given by 

S= I \m'(z)\ 2 dz = 2 [ F(m(z))dz (5.8) 



We expand / = F' as 

N 



f(mW) = f{fh) + f{fh)Y,^{K + <t>n) 

n=l 



+ 



(5.9) 



Our goal is to show that we can choose the hj and <f>j so that m^ N ^(r) is a high 
order approximate solution of <\5.1\i . At the j-th stage, to determine <pj, we will solve an equation 

in L°°(R), and to determine hj, we will solve an equation in L 2 (M). The Fredholm criterion 

for solvability of this equation will relate \ij to and the constraint equation (|5.3|) will then 
determine 

5.2 The first order approximate solution 

To see how this goes, we insert 



{1) _ ( d(x,rU) 

m y '(r) = m I 



+ \{hi + <j) 
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into ()5.1|) . collect terms by order in A, solve the resulting equations to find r±, fix, hi and 0i. 
The result is: 

-fh" + F'(fh) + A (-K^fti - ti( - <t>'{ + F"(fh)hx + F"(fh){4>x + hi) + 0(A 2 ) . 



Since, by the definition of m the term multiplied by A vanishes, we need to equate to the 
coefficient of the term of order A. This corresponds to put 

-K^fti + {-ti[ + F"(fh)hx) + (-0i + F"(m)4>x) + ni = . (5.10) 

Define the operator £ by 

£~£ + f<«». 

Also define £o by 

d 2 d 2 1 

A) = ~ + i^(i) = ~ + -. 

In the present example, 1/x = 2, and so /2™ 1 is given by convolution with the Helmholtz Green's 
function 

2^/2 

£q 1 maps constants into constants and preserve the parity properties of functions. 
Then we can rewrite (|5.10j) as 

Ch x = K^fti + (F"(l) - F"(rh))4>x - - Mi • (5.11) 

The first two terms on the right hand side are rapidly decaying in z. In order for the entire right 
hand side to rapidly decay, and thus belong to L 2 (R), we require that the last to terms cancel. 
That is, we require 

C (pi + nx = • 

This is solved by 

<Pi = -X/xi , (5-12) 

and so (|5.11|) reduces to 

Ch x = K {l) m + (1 - X F"{m))nx • (5.13) 

We have ^ 

<f>x = Co Vi = /"i (5-14) 

X 

The equation for hi becomes 

Chx = K^m' + (1- X f'(m)fix = (5.15) 

The null space of C is spanned by fh' , and so the Fredholm criterion says that 1)5. 15|) is solvable 
if and only if K^fh 1 + (1 — ;^/'(m))/ii ls orthogonal to to'. Multiply by fh' and integrate. Using 
the fact that J R f'(fh)fh'dz = 0, and J* R fh'dz = —2, we obtain, 



K {1) (J (m'fdz 



-2/ii • 
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We see that to leading order in A, the curvature must be constant, and so T^ N ' must be a circle, 
since we are considering values of n that are close to —1. Let be the radius of the circle, which 
is, as yet, undetermined. 

Using l|5.8|) to express the integral in terms of S, we obtain 

K^S , , 

/ii = — . (5.16) 

Moreover, since fh(z) = — tanh(z/\/2), it is easy to see that both fh'(z) and (1 — xfi™)) are 
proportional to sech 2 (z/\/2), and so with m given by (|5.16j) . the right hand side of 1)5. 1 5 j) vanishes 
identically, and we see that hi = 0. From (|5.14l) and (|5.16|) . we have 

K^S 

4>i = X—j- • (5-17) 

Finally, we determine using the approximate constraint 1)5. 3|) . Since = we have 

that 

(1) _ / d(x,rW) \ x s 

m y ' = m + 



A / 2r« ' 
Toward this end, note that by (|5.7j) . 

'd{x,VW)\ (rW) 2 



in* 



Thus, 



^m«(x)dx = -l + 2vr^j- + A</»i + 0(A^ 



n + 0(A 2 ) = -l + -g- + 0(A 2 ) 



(5.18) 



This yields 



2vr[(r( 1 )) 2 -l] + A|O A |^ y = 0, (5.19) 

which is a cubic equation determining rW. One could solve it explicitly - it is, after all, a depressed 
cubic. However, the important point is that the of the three roots, only two are positive, and the 
largest one is exactly the radius determined by the Biskup-Kotecky-Chayes prescription. 

To see this, we compute the free energy FijrS 1 ^). First note that, with the scaling we used, we 
have 

r A 2 



F(mF>) = A~ 2 / [— |Vm (1) | 2 + F{m (1) )]dx . 
By using the expression of and Q5.7|) to pass to the variables (s, z), we get 
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From (|5.18|) and (|1.6|) we compute 



_ 2vr(l-(rW) 2 ) 

Jl " — XinXj + 0(A) - 



Therefore 

JP( m (i)) = A" 1 ^TrrWfi + 2vr 2 (l - (r«) 2 ) 2 ) + O(l) 

This is exactly the phenomenological free energy of [I] for the Free energy functional 1)1. Also, 
the Euler-Lagrange equation for it reduces to (|5.19|) , Hence m^ 1 ', with chosen to be the 
minimizing solution of (|5.19j) is essentially exactly the trial function we used to obtain the upper 
bound in Section 2. The only difference is a slight adjustment of the additive constant so that the 
constraint in (|1.2|) is exactly satisfied. 

The crucial feature in the free energy function that is responsible for this is that in this 
case we found h\ = 0. 

5.3 The second order approximate solution 

Going on to second order is not difficult since h\ = 0, and 0i is an explicit constant. The next 
order displays some new features, but once the calculations are carried out to second order, it will 
be clear how to go on to arbitrary order. 
The equations for hi and 02 are 

-^(^2 + <h) - K^m' + f'(m)(h 2 + &)~/"(ro)(/n + X ) 2 + fi 2 = 
As before, write 

£02 = £ O 02 + (f'(fh) - /'(1))0 2 , 

and we then have 



Ch 



2 



K Wrn' + (F"(m) - /'(!))&] - A)02 - ^/"(m)0 2 - fi 2 ■ 



The terms in square brackets are rapidly decaying as long as is bounded. To eliminate the other 
terms, we choose 

A)02 = --F"'(fh)(t)l - I2 2 • 

This is easily solved in L°°(M) using the Helmholtz Green's function. With this choice of 02, we 
determine h 2 through 

Ch 2 = -K®m' - (F"{m) - F"(1))0 2 . 

The right hand side must be orthogonal to to'. Recall that (/'(to) — /'(l)) is even, and in fact 
is a multiple of to'. Since /"(to) = 6to, which is odd, only the even part of 02 is relevant in the 
solvability condition. But clearly, the even part of 02 is given by 

(02)even = ~X^2 ■ 

Thus, the solvability condition becomes 



( [ (m'(^)) 2 d^) + - X 2fi 2 = . 
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As before, with this choice, 

K^fh' + (/'(m) - / , (l))(0 2 )evcn = , 

since it is a multiple of m! , and the equation for h 2 reduces to 

Ch 2 = -{f'(m)-f'{l)){fo) odd . 

The mass condition is 

\hi m(2){x)dx = - i+2 "iSf +a ^) +a2( ^— +o(a ' 



2"7T 

= n + o( A 3) = _i + _ + o(A3). 



(5.20) 

By choosing = r^ + Xr 2 , recalling that (fa) even = we get 

xs 



T2 = ~ 



The even part of <f>^ is the most significant term among the new correction. However, it is a 
constant, so that even keeping the most significant term at second order, we still have a trial 
function of the simple type considered in Section 2. For this reason, one can expect that the upper 
bounds obtained in section 2 are in fact sharp not only in the leading order, but in the first two 
orders in powers of L. 

The procedure can be continued along the same lines to higher order thus producing approximate 
solutions to the Euler-Lagrange equations. However, what is probably more significant is that it 
can be applied to other free energy functional with non local interactions for which it will not be 
the case that hi = 0, and hence the construction of a suitable trial function is not so simple. 
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